Electronic structure theory of the hidden order material URu2Si2 
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We report a comprehensive electronic structure investigation of the paramagnetic (PM), the large 
moment antiferromagnetic (LMAF), and the hidden order (HO) phases of URu2Si2. We have 
performed relativistic full-potential calculations on the basis of the density functional theory (DFT), 
employing different exchange-correlation functionals to treat electron correlations within the open 
5/-shell of uranium. Specifically, we investigate — through a comparison between calculated and 
low-temperature experimental properties — whether the 5/ electrons are localized or delocalized in 
URu2Si2. The local spin-density approximation (LSDA) and generalized gradient approximation 
(GGA) are adopted to explore itinerant 5/ behavior, the GGA plus additional strong Coulomb 
interaction (GGA-f t/ approach) is used to approximate moderately localized 5/ states, and the 5/- 
core approximation is applied to study completely localized uranium 5/ states. We also performed 
dynamical mean field theory calculations (LDA-)-DMFT) to investigate the temperature evolution 
of the quasi-particle states at 100 K and above, unveiling a progressive opening of a quasi-particle 
gap at the chemical potential when temperature is reduced. A detailed comparison of calculated 
properties with known experimental data demonstrates that the LSDA and GGA approaches, in 
which the uranium 5/ electrons are treated as itinerant, provide an excellent explanation of the 
available low-temperature experimental data of the PM and LMAF phases. We show furthermore 
that due to a materials-specific Fermi surface instability a large, but partial, Fermi surface gapping of 
up to 750 K occurs upon antiferromagnetic symmetry breaking. The occurrence of the HO phase is 
explained through dynamical symmetry breaking induced by a mode of long-lived antiferromagnetic 
spin-fluctuations. This dynamical symmetry breaking model explains why the Fermi surface gapping 
in the HO phase is similar but smaller than that in the LMAF phase and it also explains why the HO 
and LMAF phases have the same Fermi surfaces yet different order parameters. Suitable derived 
order parameters for the HO are proposed to be the Fermi surface gap or the dynamic spin-spin 
correlation function. 

PACS numbers: 71.20.-b, 71.27.-(-a, 74.70.Tx, 74.20. Pq 



I. INTRODUCTION 

Over the past 15 years the concept of "hidden order" 
has evolved to describe the emergent behavior of various 
quantum or strongly correlated materials where the or- 
der parameter of a clear phase transition along with its 
elementary excitations remain unknown. Often modern 
microscopic measurement techniques of diffraction (neu- 
trons or photons), nuclear magnetic resonance (NMR), or 
muon spin-rotation (/iSR), etc. are unable to detect and 
characterize the new ordered phase. Yet the thermody- 
namic and transport properties unambiguously discern a 
novel state of matter appearing at a sharp transition tem- 
perature. Within this state additional unconventional 
phases may form depending upon varying parameters 
such as pressure, magnetic (electric) fields, and doping. 
Although there is at present no comprehensive review of 
the generic hidden order (HO) problem and its relation 
to quantum criticality, the HO concept is beginning to 
make headway into the recent literature.'i"Sl 

A prototype system for this behavior is the intermetal- 
lic compound URu2Si2, discovered 25 years ago.l^i^This 
material displays strong electron correlations such that 
the U 5/ magnetic moments are dissolved into hybridized 
bands near the Fermi surface and a moderately heavy 



Fermi liquid forms at temperatures below ca. 70 K !^" ! ^^ ! 
Then at 17.5 K the HO state appears via a dramatic 
(second-order) phase transition.'^'EI The above-mentioned 
techniques fail to discern the order parameter and cannot 
characterize its elementary excitations. Great attention 
has been devoted to studying this system with the aim of 
uncovering its hidden nature. A vast collection of exper- 
imental data has been gained and excellent single crys- 
tals are now available for definitive investigations thereby 
eliminating extrinsic effects of impurities and stress (see, 
e.g., Ref . fT2l) . In addition there are numerous theoretical 
proposals and exotic models spanning many yeara^^'Sll 
that have, however, not come to full grips with many 
aspects of the experimental behavior. 



Recent investigationgl^EllMl on good single crystals 
have mapped out the phase diagram of URu2Si2. Apart 
from the paramagnetic (PM) phase and the HO phase be- 
low 17.5 K at ambient pressure, there is also the large mo- 
ment antiferromagnetic (LMAF) phase, which appears 
with modest pressure of about 0.5 GPa and is charac- 
terized by uranium moments of 0.4 /i^ m a typc-I anti- 
ferromagnetic (AF) arrangement. Surprisingly, the bulk 
properties of the HO and LMAF phases are very much 
alike. Very similar, continuous changes in the thermo- 
dynamic and transport quantities have been reported 



for both phases fi^MElElM] ^ comparable Fermi surface 
gapping occurs for the transitions from the PM phase to 
the HO and LMAF phases, respectively. This similarity- 
which has been called adiabatic continuity^SL extends 
to the Fermi surfaces of the HO and LMAF phases. 
De Haas-van Alphen experiments detect no significant 
differences between the Fermi surfaces of the HO and 
LMAF phases^ and, consistently, neutron scattering ex- 
periments find the same nesting vectors.'^ Nonetheless, 
the HO and LMAF unmistakably have different order pa- 
rameters; simple magnetic order in the LMAF phase but 
an unknown order parameter in the HO phase. Neutron 
and x-ray scattering experimentd^^MS] detected a small 
magnetic moment ~0.03 iib m the HO phase, but this 
small moment is currently considered as a parasitic mo- 
ment that is not intrinsic to the HO phase. ^^ '^^ ^"^ Other 
differences between the HO and LMAF phases is that 
below 1.2 K and only out of the HO an unconventional'*^ 
superconducting state appears, which is the subject of 
recent interest .HSHm A further salient difference between 
the HO and LMAF phases is that inelastic neutron ex- 
periments detected a mode of AF spin-fluctuations in the 
HO phase which freezes to the st atic a ntiferromagnetic 
Bragg peak in the LMAF phase .i-^^ SUSS] 

As a starting point towards a full theoretical under- 
standing of the intriguing electronic structure of URu2Si2 
state-of-the-art band structure calculations are required. 
We present here detailed investigations of the electronic 
structures of the PM and LMAF phases, using various 
computational methods. On the basis of the obtained 
electronic structures, we analyze in how far the known 
physical properties of URu2Si2 can be explained from 
these underlying electronic structures, and draw conclu- 
sions on what the valid electronic structure of URu2Si2 is, 
emerging from the electronic structure calculations. Sub- 
sequently, we focus on the implications for a prospective 
explanation of the HO. Also, we expand on the "dynami- 
cal symmetry breaking" model for the HO, which we have 
recently proposed.l^ Details of this model are given and 
we relate the model to a larger collection of experimen- 
tal properties. We also compare the derived electronic 
structure and the HO model to other recent proposals. 

In the following we first consider an issue that is central 
to the current discussion of model explanations of the HO 
phase. 



II. ITINERANT OR LOCALIZED 5/ 
BEHAVIOR? 



One of the most intriguing questions regarding the 
electronic structure of URu2Si2, and consequently the ex- 
planation of the HO, is whether the uranium 5/'s are lo- 
calized or delocalized. Single-ion theories of the HO, such 
as, e.g., quadrupolar or octupola r ordering, are bas ed 
on the assumption of locahzed s/'s.lSIMIMmiMMIllHIl 
This important issue of the degree of 5/ localization has 
been controversially discussed recently. Several theories 



adopt the picture of localized 5/ states from the out- 
set, however, an examination of the grounds for this is 
needed. A thorough examination seems to unveil that 
there is little compelling experimental evidence for local- 
ized 5/'s. Smoking-gun evidence for localized / states 
would be the classical observation of crystal electrical 
field (CEF) / excitations in neutron experiments, but 
only itinerant spin excitations have been detected and 
CEF excitations have never been observed for URu2Si2 
(see, e.g., Ref . [53|) . Another indication of a CEF excita- 
tion could, e.g., come from measured specific heat curves, 
in which humps or peaks could signal the occurrence of 
CEF excitations. The measured C/T curve of URu2Si2 
shows a maximum at 70 Kp' which has sometimes been 
interpreted as evidence for a CEF transition. However, 
later measurements^ of the C/T of ThRu2Si2, which has 
no occupied 5/'s and hence no CEFs, revealed a very 
similar maximum at the very same temperature. This 
suggests that the peak at 70 K is more likely related to 
the same underlying lattice structure and not to a CEF 
transition of 5/ states. The shape of the measured mag- 
netic entropy Sm{T) in the PM state does not corres pond 
to Schottky-type anomaly expected for CEF levels.'^^Ml 
Also, very recent scanning tunneling spectroscopy (STS) 
measurements could not detect any CEF splitting of the 
5/'s.^^ Consistently, the susceptibility of URu2Si2 does 
not show Curie- Weiss behavior near the HO temperature 
that might indicate localized / states, r athe r Curie- Weiss 
behavior commences only above 150 K.f^EZl 

One particular piece of experimental evidence in favor 
of localized 5/'s has come from inelastic neutron scatter- 
ing experiments,^^ in which a small inelastic peak was 
observed at 363 meV. This peak has been interpreted as 
a signature of an inter-multiplet transition.'^ A similar 
peak has been observed for UO2, which is indeed known 
to have a localized 5/^ configuration. However, for UO2 
CEF excitations, too, were definitely observed with in- 
elastic neutron scattering (see, e.g., Ref. l59|) . Inelastic 
neutron experiments^- also detected a small peak at 363 
meV for ThRu2Si2, which indicates that the peak might 
not be due to an inter-multiplet excitation. In addition, 
a similar peak has been observecP^l for URhAl, which 
is however known to be an itinerant 5/ material.^* The 
origin of the neutron peak at 380 meV in URhAl has 
consequently been debatedpl the issue is not completely 
solved, but it could be an artifact related to the measure- 
ment apparatus. 

Several other experimental data rather advocate the 
existence of delocalized 5/ electrons in URu2Si2. High- 
resolution photoemission spectroscopy (PES) using He I 
and He II radiation gave evidence for a typical delocal- 
ized 5/ response in the He II— He I difference spectrum.'SS 
A similar difference spectrum has been observed for itin- 
erant U-metal and UGas.l^In addition, angular resolved 
photoemission spectr oscop y (ARPES) revealed disper- 
sive bands in URu2Si2,'2^'^yet it still needs to be clarified 
what the dominant character of the observed bands is (/ 
related or not). On the other hand, very recent He I 



ARPES measurements provided a picture of an almost 
flat band which sinks through Ep at the HO transition.- ^ 
The picture of a narrow band very close to Ep may how- 
ever arise from the special data treatment, i.e. division by 
the Fermi function and double-differentiation technique) 
which always tends to give an impression of a flat state 
near Ep. 

The AF phase of URu2Si2 is commonly referred to as 
the large moment antiferromagnetic phase. This name 
suggests that the 5/'s in the LMAF phase might be par- 
tially localized. However, in spite of its name, the ura- 
nium moment in the LMAF phase is actually relatively 
small and not the typical moment of a localized 5/ ma- 
terial. For example, the 5/ states of the cubic uranium 
salt USe are known to be closer to 5/ localization, but 
still exhibit some f — d hybridization, which leads typi- 
cally to spin, orbital, and total moments of -1.1, 3.1, and 
2.0 /is, respectively, for USciS^ The total moment on U 
in AF URu2Si2 is with 0.4 fiB quite far from such value. 
Instead, the U moment is much closer to values of 0.6 ub 
measured for an itinerant 5/ material such as UGas.'^ 

It also deserves to be mentioned that the results of 
recent positron annihilation experiments on URu2Si2 
proved to be incompatible with localized /'s, but are in 
contrary in good agreement with delocalized f^sP^ Also, 
recent neutron scattering experiments detected itinerant 
5/ spin excitations.^^ 

Altogether, there does not appear to be clear, com- 
pelling evidence for localization of 5/ electrons in 
URu2Si2. On the other hand, there exists a body of 
evidence in favor of delocalized 5/'s. Nonetheless, the 
decision on localized-itinerant behavior should be con- 
cluded from an extensive comparison of calculated and 
experimental properties, which will be presented below. 

In the following Section we first outline the here-to- 
be applied first-principles based techniques to study the 
electronic structure of URu2Si2. With these different ap- 
proaches we can treat the full range of 5/ behavior, from 
delocalized to localized. In view of the above considera- 
tions regarding the itinerant or localized 5/ behavior, our 
main focus will be on the 5/ band description. The ap- 
plied DFT-LSDA and LSDA+C/ approaches can provide 
only ground-state T = K properties. The tempera- 
ture dependence of quasi-particle spectra will be treated 
through dynamical mean field theory (DMFT) calcula- 
tions. 



III. COMPUTATIONAL METHODOLOGY 

Our calculations are based on the DFT as well as on the 
DMFT. Specifically, for the treatment of the DFT static 
exchange-correlation potential we have employed the lo- 
cal spin density approximation (LSDA),Eil the general- 
ized gradient approximation (GGA),l^and also orbital- 
dependent extensions (LSDA+f7, GGA+U) to include 
the influence of strong on-site Coulomb correlations. 

In our calculations we have used three accurate full- 



potential, relativistic electronic structure codes. These 
are the full-potential local orbitals (FPLO) method^^lZil 
and the full-potential linearized augmented plane wave 
(FLAPW) method; the latter we employed both in the 
WIEN2K[Z51 and Kansai implementations. We have ver- 
ified that the three codes give, on the self-consistent 
LSDA level, identical results for the electronic structure 
of URu2Si2. 

In the FLAPW calculations the relativistic valence 
states are computed within the full, non-spherical po- 
tential. The relativistic spin-orbit interaction (SO) was 
included self-consistently,'^ and, in the WIEN2K cal- 
culations, we used the relativistic local orbitals exten- 
sion of the scalar-relativistic FLAPW basis to treat 
accurately the 2pi/2 semi-core states.^^ The product 
of Rmt and maximum reciprocal space vector (Kmax), 
i.e., the basis size determining parameter [RKmax) was 
set to 7.5 and the largest reciprocal vector G in the 
charge Fourier expansion, Gmax, was equal to 12. We 
used about 5000 fc-points for self-consistent convergence. 
With the WIEN2K calculations we have employed the 
orbital-dependent GGA-I-C/ method with around mean- 
field double-counting correction,'^ in which an additional 
on-site Coulomb interaction, expressed by the Hubbard 
U and exchange J parameters, is introduced for the 5/- 
states manifold. 

In the relativistic full-potential FPLO calculations^ 
the four-component Kohn-Sham-Dirac equation, which 
implicitly contains spin-orbit coupling up to all orders, 
is solved self-consistently. We used in the FPLO calcu- 
lations the following sets of basis orbitals: 5/; QsQpQd; 
7s7p for U, 4s4p4d; 5s5p, and 3s3p3d, for Ru and Si, re- 
spectively. The high-lying 6s and 6p U semicore states, 
which might hybridize with other valence states are thus 
included in the basis. The site-centered potentials and 
densities were expanded in spherical harmonic contribu- 
tions up to Imax — 12. Brillouin zone (BZ) sampling was 
performed with maximally 20x20x20 fc-points. 

For the DMFT calculations we have used a recently 
developecP^ full-potential, relativistic LSDA-hDMFT 
method. For a detailed review of the DMFT method we 
refer to Ref. [79l In our DMFT calculations we use the 
spin-polarized T-matrix fluctuation-exchange (FLEX) 
impurity solver^*' ^^ for generating the self-energy. This 
impurity solver is expected to be applicable to moder- 
ately correlated materials, as, e.g., uranium intermetal- 
lic compounds. The local Green function is computed 
employing Kohn-Sham states which are obtained from a 
relativistic LSDA-f C/ calculation. The Coulomb U and 
exchange J quantities are connected to the two-electron 
integrals of the Coulomb interaction of the / electrons 
through the effective Slater integrals i^K (0 < k < 6), 
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J ^ 286F2-H^95f4+250F6 _ OveraU self-consistency has been 
achieved through iterative feeding of the density matrix 
of the local Green function in the next LSDA-I-C/ loop, 
and back feeding of the new solutions in the DMFT loop. 
In the DMFT calculation we used 8192 Matsubara fre- 



quency points to compute the temperature dependent 
quasi-particle spectrum, which was obtained using a Pade 
approximation to the frequency dependent lattice Green 
function G{k, iuj). Within the present DMFT FLEX im- 
plementation temperatures down to about 100 K can only 
be be reached. 

In our investigations we employ the body-centered 
tetragonal (i.e., ThCr2Si2) structure with space group 
No. 139 for paramagnetic URu2Si2 and the simple tetrag- 
onal structure with space group No. 123 for AF URu2Si2. 
The ST unit cell volume of AF URu2Si2 is twice the bct 
unit cell volume of PM URu2Si2. The space group of 
the HO phase has not yet been definitely established, be- 
cause the symmetry breaking in the HO phase is as yet 
to be unveiled. In Fig. [TJthe Brillouin zones of the bct 
and ST structures are shown with high-symmetry points 
indicated. The ST BZ corresponds to a folding of the bct 
BZ at ±iZ. In the bct BZ we have additionally labeled 
several non-high-symmetry points (S, F, Y, and A) for 
later discussion. 




FIG. 1: The Brillouin zones of the body centered tetragonal 
phase (space group No. 139, left) and the simple tetragonal 
phase (space group No. 123, right) with high-symmetry points 
indicated.'^ 



IV. RESULTS 
A. DFT delocalized 5/ electron calculations 

1. Structural optimization 

To start with, we consider the structural properties 
of URu2Si2, that is, the equilibrium lattice coordinates, 
bulk modulus, and equation of state. Several experimen- 
tal investigations of the structural properties of URu2Si2 
have been reported.'^'^SHlSl Por comparison to the avail- 
able data, we have performed ab initio optimizations of 
the equilibrium volume, the c/a ratio, and the internal 
Si coordinate, zsi. These optimizations have been per- 
formed on the LSDA level, both for the PM and LMAF 
phases. In Fig. |2] we show the computed total energy 
versus unit cell volume. Both PM and LMAF total 
energies are given for the double unit cell, to convene 
comparing the two phases. The theoretically predicted 
equilibrium unit cell volume is about 1.7% smaller than 
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FIG. 2: (Color online) The ah initio computed total energy 
versus volume for the PM and LMAF phases of URu2Si2 . The 
vertical dashed line indicates the experimental equilibrium 
volume,'^ the inset shows the ah initio computed pressure 
dependence of the unit cell volume for the PM and LMAF 
phases, normalized to the zero-pressure volume Vo, together 
with experimental data points of Ref. 1851 



the experimental volumes, which are 162.9 A'^ (Refs. [7] 
and "531), respectively, 162.6 A^ (Ref. |MD- Hence, the 
theoretical value is in very good agreement with experi- 
ment. As Fig. [2] illustrates the total energies of the PM 
and LMAF phases are very near one another. The to- 
tal energy of the LMAF phase is computed to be only 7 
K per formula unit deeper than that of the PM phase. 
This is in itself a remarkable finding, which appears to 
be a specific feature of URu2Si2. The inset of Fig. [2] 
presents the computed volume versus pressure depen- 
dence of URu2Si2. With pressure the antiferromagnetic 
state becomes slightly more stable. The inset includes re- 
cent experimental data points of Ref. [85l A fit of the com- 
puted volume versus pressure curves gives a bulk modu- 
lus Bo of 204 GPa (208 GPa) for the LMAF (PM) phase. 
The recent pressure experimenlpSI obtained a value of 190 
GPa; an older experiment reported a value of 230 GPa.^ 
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FIG. 3: (Color online) Total energy optimization of the cja 
ratio of URu2Si2 for the PM and LMAF ph ases . The dashed 
vertical line denotes the experimental value! ^ * ^^ ^ 

The optimized theoretical cja ratio is shown in Fig. p] 



The obtained theoretical cja ratio almost coincides with 
the experimental value (2.32)!2'S3lThe optimized cja ra- 
tio of the LMAF phase is found to be just a small fraction 
larger than that of the PM phase. We note that x-ray 
diffraction experiments have been unable to detect any 
difference i n the lattice constants of the PM, LMAF, and 
HO phasesl^^EZI Only dilatation experiment s'^ ^'^'^' could 
so far detect tiny differences in both the a and c lattice 
constants of the three phases; the c axis lattice constant 
of the LMAF and HO phases are elongated with a few 
parts in 10~^, as compared to the PM phase (at higher 
temperature). The a axis of the HO and LMAF phases is 
contracted by a few parts in 10~^. As a result, the cja ra- 
tio increases^ with about 10^"' from the PM above 17.5 
K down to the LMAF phase at 10 K. The optimized cja 
ratio of the LMAF phase is consistently computed here 
to be about 10^^ larger than that of the PM phase, in 
agreement with experiment.'^ We refrain however from 
a more detailed comparison because we cannot make a 
meaningful quantitative statement for such tiny numbers. 
Using the optimized c axis lattice constant and vol- 
ume, the theoretical a axis lattice constant is about 0.6% 
smaller than the experimental value .I^EMSl As the LSDA 
approach is know to produce a small overbinding, the cor- 
respondence with the experimental lattice constant can 
be regarded as very good. 
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FIG. 4: (Color online) Total energy optimization of the special 
Si z-position in URu2Si2. The zero position of zgi is taken at 
0.371 (Ref.,83,). 

The BCT structure of URu2Si2 has one internal coordi- 
nate, the Si z-position. Results for the total-energy op- 
timization of the zsi coordinate are given in Fig. |4J The 
theoretical LSDA value is found to be 3% smaller than 
the experimental value of 0.371 (Ref. [55)) . A more re- 
cent experimenlF^ obtained a somewhat smaller zgi value, 
0.3609, which would agree quite well with the theoreti- 
cal result. As will be shown in more detail below, the 
essential physical properties of URu2Si2 are stable with 
respect to moderate variations of the unit cell dimensions 
and the internal zsi coordinate. 

Altogether, the ah initio structural optimization shows 
that the crystallographic properties of URu2Si2 are well 



described by the LSDA approach, which intrinsically is 
based on the assumption of itinerant 5/ electrons. It is 
known from computational investigations for other ac- 
tinide materials that, when these have 5/ states that are 
localized, the LSDA approach usually does not provide a 
good description of the lattice properties (see, e.g., Ref. 
[55|) . Such deviant behavior is not found here for URu2Si2. 



2. Energy band dispersions 




FIG. 5: (Color online) Computed LSDA energy dispersions 
of URu2Si2 in the PM and LMAF phases, both are shown, 
for sake of comparison, in the simple tetragonal BZ. 

We first consider the outcome of LSDA 5/-itinerant 
calculations for the electronic structure of URu2Si2. In 
Fig. [5] we show the computed LSDA energy dispersions in 
the PM and LMAF phases for the experimental lattice 
parameters;^-' To draw a comparison, both sets of dis- 
persions are given for the double unit cell (space group 
No. 123). As has been noted recently by us, the en- 
ergy dispersions of these two phases are very similar.l^S 
The dispersions of the AF phase are almost on top of 
those of the PM phase, except for some influence of the 
exchange splitting of 5/ related bands. This finding cor- 
roborates fully with the compute tiny total energy differ- 
ence between these two phases. A degeneracy of cross- 
ing bands occurs near the Fermi level, as can be rec- 
ognized along the F— M and X— F symmetry directions. 
The degenerate band crossing, existing in the PM phase 
along the F— M direction just below Ep, is lifted in the 
LMAF phase, due to a re-hybridization, and thereby a 
small gap opens. A similar degenerate crossing point 
along the X— F direction is however not removed in the 
LMAF phase. Through a larger part of the BZ degen- 
erate crossings of the two bands exist, yet the opening 
of a gap in the AF phase does not happen uniformly 
over the Fermi surface (FS).'22 This gapping is related 
to a FS instability of URu2Si2 in the PM phase, where 
degenerate band crossing ( "Dirac points" ) occur off the 
high-symmetry directions, between the F— M and F— X 
directions. These degenerate points are removed in a 
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FIG. 6: (Color online) Enlarged view of the PM and the 
LMAF bands along the F— M (E) high-symmetry direction in 
the simple tetragonal BZ. The character of the LSDA bands 
is given through the color of the bands (green: Ru 4d, red: U 
5/ in LMAF phase, blue: U 5/ in PM phase). The amount 
of Ru 4d or U 5/ character in the respective bands is given 
by the thickness of the bands. 



transition to the LMAF phase, leading to a fc-dependent 
FS gapping that is largest in the z — plane.I^ 

An enlarged view of the PM and LMAF energy bands 
along the F— M direction is shown in Fig. |6] In addition 
we have highlighted the orbital character of the bands 
through the colors and the amount of orbital character 
through the thickness of the bands. Ru Ad character is 
shown by the green color; the bands that consist primar- 
ily of Ru Ad character appear about 0.20 cV below Ep. 
The bands closer to the Fermi energy contain dominantly 
U 5/ character, as is shown by the blue color in the PM 
phase and the red color in the LMAF phase. A small 
admixture of Ru d character is nevertheless present. The 
lifting of the degenerate band crossing is clearly borne out 
uranium 5/ dominated states. The bands with mainly 
Ru Ad character are unaffected. The gap opening due 
to the re-hybridization of states in the LMAF phase is 
about 60 meV wide and located in a narrow reciprocal 
space region at a distance of about 0.25a* — 0.3a* from 
the r point. These values agree well with those observed 
in recent STS measurements.'^ 

URu2Si2 is known to be a compensated metal (see, 
e.g., Refs. \M and [1D|. The LSDA-computed energy 
bands of URu2Si2 (Fig. [s]) are fully consiste nt wit h the 
property, as has been pointed out recently r'^ * ^° l Both 
the opening of the gap in the symmetry-broken phase 
and the compensated metal character are closely con- 
nected to the uranium 5/ occupancy. Our LSDA calcu- 
lations predict a 5/ occupancy of 2.7,^^ a value which is 
consistent with recent electron energy loss-spectroscopy 
(EELS) measurements. fSSl 

The lifting of a FS instability in the LMAF phase is 
a significant feature of URu2Si2 obtained from ab initio 
calculations. The gap appearing around the Fermi level 
is narrow and might therefore sensitively depend on the 
lattice constants. In Fig. [7] we show the influence of the 




FIG. 7: (Color online) Computed LSDA energy dispersions 
of antiferromagnetic URu2Si2 for several unit cell volumes 
around the experimental volume (162.9 A^). 



volume on this feature. For a range of volumes about the 
experimental volume the gapping property is found to 
be stable. We have similarly investigated the influence of 
the zgi coordinate on the FS gapping (not shown). Also 
for the zgi coordinate we find that the gapping property 
is stable for a range of values around the experimental 
one. 

To end this LSDA/GGA band structure section we 
briefly mention that several LSDA electro nic st ructure 
calculations have been reported for URu2 Si2 j ^^ l ^'^I Eil Roz- 
ing et al^ and Ohkuni et alW^ reported LDA calcu- 
lations for PM URu2Si2, Yamagami and Hamada^Sl re- 
ported LSDA calculations of antiferromagnetic URu2Si2. 
The non full-potential calculations^" ^* for the PM phase 
are in reasonable agreement with our full-potential re- 
sults. Our energy bands and FS of AF URu2Si2 are how- 
ever distinctly different from earlier published results.!^ 
A reason for this difference is not known. As mentioned 
before, we have verified that independent state-of-the-art 
electronic structure codes give nearly identical results. In 
Ref. [SSJan AF state with nearly compensating antiparal- 
lel spin and orbital moments (possibly obtained with an 
orbital polarization term) is proposed as a solution for 
the small moment antiferromagnetic (SMAF) phase; but 
the SMAF phase is nowadays considered to be parasitic 
rather than intrinsic.SS 



3. Density of states 

The computed total and partial densities of states 
(DOSs) of URu2Si2 are shown in Fig. [SJ Note that the 
Fermi level (at eV) falls precisely in a sharp minimum 
of the total DOS. The contribution of the uranium 5/ 
states increases, starting from about 2 eV below Ep; at 
Ep and up to 1.5 eV above Ep, the uranium 5/ dominate 
the DOS. The hybridized Ru Ad states extend from -6 eV 
to 3 eV; the hybridized Si pstates extend over the same 
energy interval. Yang et a/.l^ performed a photoemission 
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FIG. 8: (Color online) Partial density of states (DOS) of 
URu2Si2 in the paramagnetic and large moment antiferro- 
magnetic phases. The inset shows an enlargement of the PM 
and AF DOS around the Fermi energy (at eV), illustrating 
the partial gapping occurring at the Fermi energy. 



study of URu2Si2, in which they employed both He I and 
He H radiation to locate the energy position of the 5/'s 
relative to the Fermi level. From the difference of the 
He H and He I emission spectra Yang et al. inferred that 
the U 5/ states are relatively delocalized and energeti- 
cally extend from 1.5 eV binding energy up to the Fermi 
energy, where they assume a maximal contribution. The 
Ru Ad states were found to be located at about 2 eV 
binding energy. The computed LSDA DOS is in good 
agreement with these findings. Fig. [8] furthermore illus- 
trates that the total DOSs of the PM and LMAF phases 
are very similar, as expected. The difference between the 
PM and LMAF DOSs is largest in the 5/ energy inter- 
val; the inset of Fig. [8] shows the total DOS of the two 
phases on an enlarged energy scale close to Ep. Due to 
the gap opening on a part of the FS in the LMAF phase, 
the DOS minimum at Ep deepens. 



4- Magnetic moments 

The total magnetic moment on the uranium atoms in 
the LMAF phase is reported to be 0.40 fj,B in recent 
measurements. 03] Our ab initio calculations give a total 
moment of 0.39 /ub, in good agreement. The respec- 
tive spin and orbital moments are Ms = 0.36 hb and 
Mi = —0.75 fiB: i-G-j the orbital moment is antipar- 
allel to the spin moment and twice as large. Detailed 
measurements of the separate spin and orbital moments 
have not been reported. Nonetheless, it can be inferred 
that the computed spin and orbital moments separately 
are in good agreement with experiment. Recent neu- 
tron form factor measurements^ indicated a value for 
C2 = Ml/{Ml + Ms) of about 1.8±0.2. From the theo- 
retical values of M/, and Ms we obtain C2 = 1.9, im- 
plying that the predicted values for the spin and or- 



bital moments are indeed consistent with experiment 
(which, re-calculating from C2 and total moment, gives 
Ms = 0.32 ± 0.04 hb, Ml = -0.72 ± 0.08 hb)- 



5. Transport properties 

The thermal and charge transport prop- 
erties of URu2Si2 are experimentall y well 
documented." 9'i"'3 0|3i | 34 | 46 | 97H Iol The normaPEESHMl 

and Hal l''^° ^^ ^" ^"'^' resistivity as well as the thermal 
conductivity and Nernst effedpSESEool display a clear 
signature of the HO transition, consistent with the 
sudden removal of a part of the FS at Tq. This distinct 
jump in the transport quantities is present both at the 
phase transition from the PM to the HO phase and 
that of the PM to the LMAF phase, l^niMIinil yet detailed 
charge transport measurements revealed that the FS 
gapping in th e LMAF phase is distinctly larger than in 
the HO phase.l^SIllIlIIioTIiol] yi^^i^ ^t al^ expressed the 

FS removal in terms of the opening of a partial FS gap 
A. The gap opening in the HO phase was measured 
to be about 70 - 80% of th at of the L MAF phase 
(A,„ « 75 K, A,„,, « 100 K).[2MIESIIimi 

The resistivity change in the transition from the PM 
to the LMAF phase is accessible from the electronic 
structures. To compute the electrical conductivities in 
these two phases, we have used the Kubo linear-response 
formulation in constant relaxation time approximation. 
Apart from the Fermi velocities, the conductivity ex- 
pression contains an unknown electron lifetime which en- 
ters as a constant pre-factor. The electron relaxation 
time dependence drops out when the resistivity change 
is evaluated. For the DFT-GGA electronic structure we 
compute an unexpectedly large and also anisotropic re- 
sistivity change due to the opening of the gap at the 
PM to LMAF phase transition. The computed resistivity 
jumps are (Plmaf-Ppm)/Ppm('^I|c) = 620%, and(ft„ - 
/9p„)//9p„(J||a) = 160%. IntheexperimentPnSMoSltije 

resistivity signal is superimposed on a large background, 
~ Po+AT"^, due to incoherent and phonon scattering. We 
have subtracted this background to obtain the resistivity 
change due to the partial FS gapping only. In this way 
we obtain the measured resistivity changes in the PM to 
HO phase transition, which are about 400% and 100% 
for current along the c, respectively, a axis.^*^- These val- 
ues are consistent with the resistivity jumps computed 
for the LMAF phase, but they are somewhat smaller. 
This might be related to the fact that the measured re- 
sistivity jump pertains to the HO phase, in which the 
partial FS ga p is sma ller, about 70 — 80% of that of the 
LMAF phase.^ '^^ * ^°^ l Hence, the estimated experimental 
resistivity changes in the PM to LMAF transition would 
be higher (a plain scaling would give 500% and 125% 
for J\\c and J\\a, respectively). We also mention that 
the conductivity calculations pertain to the T = K 
coherent electronic structure, whereas the measurements 
were performed in the temperature range around To. The 



observed anisotropy ratio of the resistivity jump is 4 : 1 
(Ref . llOSp , a value which is in very good agreement with 
the theoretical anisotropy ratio of 3.9 : 1 predicted on the 
basis of the DFT-GGA electronic structures. 

Magneto-transport studies revealed that UR u2Si2 is 
a low-carrier, elect ron-ho le compensated metal.HSlSniinni 
As was pointed outP^Eni recently, itinerant 5/ calcula- 
tions (LDA or GGA) indeed accurately predict this fea- 
ture for URu2Si2. The electron and hole Fermi vol- 
umes in the PM phase cancel each other within 2%.E2l 
The number of holes has been determined from Hall ef- 
fects measurements to be 0.017< nh< 0.021 per U-atom 
in the H O phas e, and 0.1 per U, respectively, in the 
PM phase .HSiSil ^g have used the computed intraband 

plasma frequency, ut^ — -;^^ — n, to determine the num- 
ber of carriers, n, in both the PM and LMAF phases 
{Vijc is the unit cell volume). In contrast to the Fermi 
volume, the plasma frequency is a FS integral and there- 
fore it counts only the carriers that contribute to the 
transport (at T = K) . The computed number of holes 
is 0.08/U-atom and 0.0185/U-atom in the PM phase, re- 
spectively, LMAF phase, in reasonably agreement with 
the experimental data. The calculated values empha- 
size that the FS gapping in the PM to LMAF transition 
strongly reduces the number of carriers by a factor of 
four. Hall effect measurements give that there are about 
five t imes less carriers in the HO phase than the PM 
phase.ESa 

The computed 5/ itinerant (LSDA or GGA) FS gap 
has recently been compared to experimental values.^* 
The FS gap at the transition to the HO phase was first 
determined by Maple et alW from specific heat measure- 
ments. This FS gap, averaged over the whole BZ, was 
estimated to be Aj^^ « 11 meV. Somewhat smaller gaps 
for the HO phase of ab out 7 meV w ere obtained from 
transport measurements, '^SEHMSIIloi] and a larg er gap of 
about 10 meV was measured for the LMAF phase.'2212^1i2il 
The FS gap which is predicted by DFT delocalized 5/ 
calculations^ is strongly fc-dependent (see Fig. p^ and 
Fig. [l6] below). The LMAF gap A^maf varies from max- 
imally 65 meV along the F— M (E) -direction to meV 
along the F— X (A) -direction. The larger theoretical 
gap obtained in certain places in the BZ is not incon- 
sistent with the smaller BZ-averaged gaps obtained from 
transport measurements. Moreover, the computed gap 
pertains to the coherent (T = K) electronic structure, 
but in the experimental analysis of the transport data 
both a k and temperature independent gap is assumed 
constant from T = K to Tq .I^EM^gMOiiios, 



6. Optical spectra 

The optical spectra of URu2 Si2 have been measured 
in several invest igat ions .'^^Stiosi xhe gapping occurrin g in 
the HO phase was observed originally by Bonn et al.^^ 
more detailed infrared optical measurements of the gap- 
ping with Re doping were performed by Thieme et al. 
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FIG. 9: (Color online) Calculated optical conductivity spec- 
tra. Re [(Taa(i^)] and Re [(Tcc(i^)], of URu2Si2. Computed spec- 
tra (at r = K) are given for different uranium total mag- 
netic moments, starting from that of the PM phase (0.0 fis) 
up to that of the LMAF phase (0.39 /^s). The inset shows 
the computed spectra on a larger energy interval. Note that 



the Re [ucc] spectra have been shifted upward by by 3 ■ 
s^^ for sake of visibility. 
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and Degiorgi et al. (Ref. il09|k the reported 
spectr#^^^nI291 are in good agre ement with one another. 



A very recent optical study^^^ obtained, however, a dif- 
ference in the low frequency response. 

Using linear-response theory we have computed the op- 
tical conductivity of URu2Si2 for the two possible geome- 
tries, E\\c and E\\a, where E is the electric field vector of 
the light. In Fig. [9] we show the calculated conductivity 
spectra. Re [o'aa(w)] and Re [crcc('^)]- The plotted optical 
conductivity spectra include both interband and intra- 
band contributions and have been calculated for a range 
of static AF moments, going from the PM phase (0.0 /is) 
up to the fuU LMAF phase (0.39 ^s). The theoretical 
spectra illustrate the effect of the increased opening of 
the FS gap in the vicinity of the Fermi energy. The opti- 
cal conductivities, being maximal for the PM phase near 
zero energy, become progressively reduced for the various 
static anitiferromagnetic phases, in particular for small 
photon energies well below 50 meV. The largest drop of 
Re [cr(w)] is obtained for the LMAF phase, where the FS 
gapping is the largest. The computed behavior agrees 
re ason ably well with experimental observations. Bonn et 
a^.'iSH measured the optical response of URu2Si2 in the 
basal plane (i?||a). They observed a reduction of the re- 
flectivity in the HO phase for photon energies below 30 
meV. Our calculation predicts a drop in Re [ac,^ below 
40 meV. Bonn et al. did not measure <Tcc {E\\c), but our 
calculations predict that a larger reduction should oc- 
cur for Re [acc\ at small energies well below 50 meV. We 
also note that together with the progressive FS gapping, 
there is a transfer of spectral weight to higher energies. 
Re [(Taa] increases slightly above 50 meV. A larger spec- 
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tral weight transfer occurs for Re [ace] for energies of 50 
meV up to 600 meV (see inset in Fig. |9]). Above 200 
meV, respectively, 600 meV, the influence of the gapping 
on the optical conductivity spectra for E\\a, respectively, 
E\\c, has vanished. 

The inset in Fig.lolshows the computed spectra for E\\a 
and E\\c on a wider energy scale. An interband peak is 
presented just above 2 eV in both Re [aaa] and Re [ace]- 
Experiment also detected a peak a t this energy.^i^ 

The experimental spectrgPSZEoll reveal a particular fea- 
ture which is not present in the calculated spectra. The 
reduction of the Drude weight at low frequencies leads 
to an increased spectral weight at 7 — 8 meV.^iS^ -^.os q^j^g 
origin of this transferred spectral weight is c urre ntly un- 
known; it was not observed in a recent study.*^^ It might 
nonetheless signal a difference between the experimental 
and computed theoretical spectra. 



7. Specific heat and magnetic entropy 

The linear-temperature specific heat coefficient of 
URu2Si2 in the HO phase is, with about 50 
m J/mol K^ pMili] not particular high, implying that 
URu2Si2 in this phase is not a heavy- fermion material. 
The S omm erfeld coefficient is comparable to t hat of, e.g., 
UGa3,C33i which is an itinerant antiferromagnet P^ ^ ^ ^^1 The 
unrenormalized specific heat coefficient calculated with 
the LSDA approach is about 9 mJ/molK^, i.e., there 
is an expected mass renormalization of six, a value not 
unusual for actinides. As a consequence, the computed 
LSDA bands will become renormalized, but not strongly. 
Our LDA+DMFT calculations (to be presented below) 
indicate a further influence of the dynamic part of the 
electronic self-energy S(a;), through which a renormaliza- 
tion of the bare LSDA band masses would occur. How- 
ever, as we can currently not compute Ile[dT,{io)/dio] 
down to low enough temperatures, we refrain from giv- 
ing values for the estimated mass renormalization. Also, 
low-energy spin- fluctuations, which are not accounted for 
in the bare specific heat coefficient, can be expected to 
give a considerable enhancement.'^ 

The entropy of URu2Si2 has drawn attention from the 
beginning.* ^ ^^- The phase transition to the HO state 
was originally discovered from a A-type anomaly in the 
specific heatpf^ the related magnetic entropy change 
in th e A an omaly is, with about 0.16 i? In 2, relatively 
large.'^EEIil Such entropy removal can, in particular, not 
be explainecJii^ by assuming a phase transition to a small 
moment antiferromagnetic (SMAF) state that at first was 
thought to be connected to the HO transition.'22Hl2l 

The total magnetic entropy Sm has been determined 
by van Dijk et aZ.E3 and Janik,E31 through subtract- 
ing the measured specific heat of ThRu2Si2, which 
has no 5/ electrons, from that of URu2Si2. From 
the specific heat difference a total electronic entropy 
SmiT) = J^{AC/T')dT' approaching i?ln4 mJ/molK 
was obtained.ISS This value is not inconsistent with our 



LSDA calculations, predicting a low-temperature 5/ 
count of 2.7. Assuming at higher temperatures an oc- 
cupancy of three 5/ electrons, a spin entropy of i?ln4 
(i.e., i?ln (25* + 1), with S = 3x 1/2) follows. 

For the HO phase, the total electronic entropy at 
To amounts to about 0.25 i? In 2 (Ref. Illip . As was 
pointed out several times, assuming the opening of a 
gap A in the electronic spectrum, the electronic specific 
heat would scale as Cm{T) ex exp(— A/fc^T), which fits 
the measur ed specific heat in the HO phase extremely 
.yyg}l| 9 | 55 |iTT] jj^ £g^g^^ ^j^g opening of gap in the magnetic 

excitation spectrum can wholly explain the entropy re- 
moved at the HO transitioni'^'^ ^^ The magnitude of the 
gap was estimated to b e about 11 meV from specific 
heat measurementspj^^ or smaller from inelastic neu- 
tron nieasurenients.--®-- 1"- Our energy band calculations 
also show that the HO gapping removes a considerable 
amount of accessible states at Ep. The bandstructure 
gap computed here is in fact larger, maximally 60 meV 
(for the LMAF phase), but it is strongly fc-dependent. 
The fc-averaged gap would thus be considerably smaller 
and be consistent with the entropy loss associated with 
the HO transition. 



B. LSDA+f/ and 5/-core calculations 




FIG. 10: (Color online) Energy dispersions of URu2Si2 in the 
LMAF phase, computed with the around-mean field GGA+U 
approach, with [7=1.4 eV and J=0.68 eV. The bands crossing 
the Fermi level are high-lighted. 

Actinidc materials with enhanced Coulomb correla- 
tions between the 5/ electrons can be computationally 
treated with LSDA-I-C/ or GGA+U calculations (see, e.g., 
Refs. 11161 11171 and IllSp , which are expected to give a 
good description for materials with a moderate degree of 
5/ localization. Actinide or lanthanide materials with lo- 
calized / electrons are conversely well described by open 
/-core calculations, in which the /'s are treated as unhy- 
bridized core electrons (see, e.g., Ref. UlSj) . 

The energy bands of AF URu2Si2 computed with the 
"around-mean field" GGA+U approach are shown in Fig. 
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10 For the Coulomb U and exchange J parameters we 
have chosen the values U = 1.4 eV and J = 0.68 eV. The 
U value can be considered as relatively small and has 
been chosen such in order not to depart much from the 
LSDA solution. Nonetheless, the computed bands in Fig. 
[To] reveal that the bands near the Fermi level are modi- 
fied to a considerable extent so that also the FS becomes 
quite different. The bands near the M point are pushed 
down, whereby new FS sheets appear. Bands near the X 
and r points are pushed upwards, whereby also a new FS 
sheet appears around X. Furthermore, two new electron 
pockets appear around A. The gap features along the 
F— M and F— X directions are strongly affected; the gap- 
ping occurring along F— M has practically vanished. As 
we shall see below, experiments support in fact the Fermi 
surface predicted by LSDA calculations. This illustrates 
that the Fermi surface and its gapping is rather sensitive 
to the Coulomb U in GGA+U calculations. This is un- 
derstandable, as the FS gap is quite small (several tens 
of meV) and the opening of the FS gap is due to a subtle 
hybridization change of 5/ bands just above and below 
the Fermi level. The Coulomb U acting on the 5/ states 
changes the 5/ band dispersions substantially. 




N L 



XY P X 



FIG. 11: (Color online) Energy dispersions of non-magnetic 
URu2Si2 computed with a localized 5/^ configuration. The 
high-lighted bands are the ones crossing the Fermi level. The 
used high-symmetry points in the bct BZ are indicated in 
Fig-0 

As mentioned before, a large number of theo- 
ries for the HO of URu2Si2 are based on the 
assumpti on of completely or nearly localized 5/ 
electrons.'HIM: i9 | 2i | 23 | 26 | 28 | 5om i| Particularly, an under- 
lying l ocalized 5/^ c onfiguration has been discussed 
recentlytSMliaililllll j^^ -^gg^f^ ^ localized 5f config- 
uration possesses very interesting properties, as is can 
sustain both a non-magnetic spin-singlet and a magnetic 
triplet configuration, something which might be related 
to the occurrence of two different phases. 



In Fig. 11 we show the energy dispersion computed for 
paramagnetic URu2Si2 with the open-core approach for a 
localized 5/^ configuration. As expected, the /-core en- 
ergy bands are very different from the bands obtained for 



PM URu2Si2 assuming itinerant 5/ valence states. Such 
band structure of URu2Si2, computed with WIEN2k in 
the PM phase was reported already in Ref. [70] and is 
therefore not repeated here. The / core energy disper- 
sions are indeed so different from the 5/ delocalized ones, 
that it makes no sense to compare them. As mentioned 
before, the 5/ occupancy obtained from LSDA itinerant 
5/ calculations is about 2.7. Even with a small depen- 
dence on the used muffin-tin sphere radius, this occupa- 
tion number is not two. Hence, it is understandable that 
very distinct energy dispersions emerge. The concomi- 
tant FS's are consequently also very different, as will be 
exemplified below when we discuss the FS of URu2Si2 in 
detail. 

LSDA-5/^ core calculations were recently also per- 
formed by Haule and Kotliar,^^ who only show a small 
reciprocal space section of the bands in a narrow energy 
interval near the Fermi level and around the F-point, yet 
their results agree with our full-potential results. In par- 
ticular, there is one band at the F point near Ep that 
has an inverted parabolic shape, and there are two bands 
with a steep dispersion crossing the Fermi level between 
F and E. 

Whether or not a localized 5/^ picture is more appro- 
priate for URu2Si2 has to be considered in the hght of 
all available experimental data. The delocalized 5/ pic- 
ture provides a quite accurate description of the known 
experimental data; this cannot be said of the localized 
5/^ configuration. In first instance one might think that 
the LMAF phase might be related to a magnetic, local- 
ized 5/ configuration, but the delocalized 5/ approach 
is thus far the only one that has provided an accurate 
explanation of the LMAF phase, which is a conventional 
antiferromagnetically ordered state without mysterious 
properties. 



C. DMFT calculations 

In our LDA+DMFT calculations we used a large num- 
ber of Matsubara frequency points (8192) when taking 
the sum over the frequencies on the temperature axis, 
nonetheless, we can only compute a finite number of fre- 
quency points and therefore the calculations are valid 
for moderately high temperatures in practice (100 K and 
above). This implies that we can investigate the influ- 
ence of dynamical electron configuration fluctuations in 
the paramagnetic phase. At high temperatures the ura- 
nium 5/ moments are expected to behave as incoherent, 
local moments. Note that the single-ion Kondo tem- 
perature is estimated to be 370 K in URu2Si2.^" With 
reducing temperature, lattice coherence between the / 
moments develops below 100 K, leading to a coherence 
temperature T* of about 70 K, which is wit nesse d by a 
maximum in the normal and Hall resistivity'^^'^ Below 
the coherence temperature T* the 5/ local magnetic mo- 
ments are incorporated into the conduction electron sea, 
which greatly enhances the electron effective masses and. 
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FIG. 13: (Color online) The LDA+DMFT quasi-particle 
bands of URu2Si2 at r=100 K (for details, see text). The 
colors indicate the magnitude of the fe-dependent spectral 
function, -:|lni G{k,E). Black lines show the LDA energy 
bands for comparison. Note that Z' denotes the Z-point in 
the neighboring bct Brillouin zone. 



FIG. 12: (Color online) The fc-integrated quasi-particle DOS 
calculated with the LDA+DMFT approach for URu2Si2 in 
the high-temperature, non-magnetic phase. 



for conventional Kondo lattice materials, is expected to 
enlarge the Fermi surface, too. 

In the LDA-I-DMFT calculations we started from the 
LDA+U approach to compute Kohn-Sham states that 
are subsequently used in the DMFT self-consistency loop. 
In the LDA-I-C7 part we assumed effective U values of 0.4 
eV and 0.6 eV (in both cases, J was set to 0.0 eV). These 
U values are chosen to approximate the more localized 
behavior of the 5/'s that is anticipated at higher temper- 
atures. 



In Fig. 12 we show the computed quasi-particle den- 
sity of states of URu2Si2 for several temperatures. Pro- 
nounced changes in the quasi-particle DOS occur around 
the chemical potential (at eV). Lowering of the tem- 
perature and increase of electron coherence leads to the 
typical opening of a quasi-particle coherence gap (also 
called hybridization gap) of about 100 meV. Concomi- 
tant with the opening of the quasi-particle gap, there is 
a build-up of spectral weight on both sides of the gap. 
The development of coherence gaps has been observed 
with infr ared optical spectroscopy for several /-electron 
materials ,^1^ but for URu2Si2 this property has not yet 
been reported. 

In Fig. [13] we show the calculated quasi-particle bands 
of URu2Si2 at T = 100 K. The bright colors depict high 
intensity of the spectral function. For comparison, the 
non-magnetic LDA bands are shown by the black lines. 
Note that the Z' point on the reciprocal space abscissa 
is positioned in the neighboring BZ. We observe that the 
LDA-I-DMFT quasi-particle bands are relatively close to 
the LDA bands. Their similarity is even more so for en- 
ergy bands below -1 eV (not shown here), because these 
bands posses less uranium / character. Some differences 
between the LDA and quasi-particle bands can nonethe- 
less be seen from Fig. [13) In the N— P— X panel the quasi- 
particle band just above Ep moves distinctly closer to Ep 
and becomes flatter. Near the X point the quasi-particle 



band below the Fermi level moves slightly upwards and 
disperses stronger downwards towards the T point. The 
fc-dependent quasi-particle DOS gives the impression of 
a band dispersing downwards from above the X point 
towards the F point and crossing Ep between the two 
points. 

Another DMFT calculation for URu2Si2 has been re- 
ported recently by Haule and Kotliar.^^ We note that 
our DMFT results are distinctly different from those of 
Ref. [5S1 We have performed self-consistent DMFT cal- 
culations using the spin-polarized FLEX impurity solver, 
starting from LDA+C/ results, which should be valid for 
the weakly correlated uranium / electrons at higher tem- 
peratures. The DMFT calculations of Ref. [511 on the 
other hand, used the one-crossing approximation solver 
together with a nearly localized uranium 5/^ configura- 
tion. The difference can be understood to arise from the 
5/ configuration used in the LDA+U band structure part 
of the LDA-I-DMFT calculation, which obviously deter- 
mines largely the gross electronic structure, whereas the 
DMFT part involving dynamical self-energy fluctuations 
induces mainly modifications of energy dispersions in the 
vicinity of the chemical potential. 

DMFT fc-dependent spectral functions can be com- 
pared to ARPES data. Several ARPES measurements 
on URu2Si2 have been reported.l^SHsZl ito et alW^ and 
Santander-Syro et alW^ both used a He I light source, 
whereas Denlinger et al)^ used tunable synchrotron ra- 
diation. Only the recent experiment of Santander-Syro et 
al. measured quasi-particle bands below Tq. These lat- 
est measurements indicate the existence of a narrow band 
just below Ep in the HO phase, as well as of an inverted 
parabolic band at fc = below Ep; the latter band was 
attributed to a surface state.'^In our 5/- itinerant LSDA 
calculation there is no such inverted parabolic band at 
the F point, but LDA localized 5/^ electron calculations 
do predict such a band (see Fig. 11), as do also the re- 
cent LDA-^DMFT calculations of Haule and Kotliar.l^ 
As photoemission at this energy is very surface sensitive, 
it could thus be that this band stems from a 5/-localized 
response of uranium atoms on the surface. Further in- 
vestigations are therefore needed to definitely establish 
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FIG. 14: (Color online) Left: LDA+DMFT quasi-particle 
bands for the directions A— N and A— P in the bct BZ 
(A =Z/2), and right: the LDA energy bands for the same 
directions. 



the origin of the inverse parabolic band. A bulk, flat 
band just below Ep at the T point is not predicted by 
our delocalized 5/ LDA or LDA+DMFT calculations. 
We note, however, first, that the LDA and LDA+DMFT 
bands have a significant dispersion along the k^ direc- 
tion, and second, that the k^ position in the BZ will, 
in an normal-emission ARPES experiment, depend on 
the energy of the used radiation. With He 1 radiation, 
a kz position between F and Z in the bct BZ will be 
probed,'^ probably being closer to F than Z. In Fig. 14 
we present computed LDA and LDA+DMFT bands for 
the midpoint, A =Z/2, between F and Z, and going to 
the P, respectively, N, high-symmetry points in the bct 
BZ (see Fig. IT]). The plotted bands illustrate that a flat 
band exists in the P— A— N plane, just below Ep, being 
mostly flat near kx, ky — 0. At his point it is still too 
early to decide whether this computed flat band does or 
does not correspond to the observed ARPES structure."^ 
High-resolution ARPES measurements with tunable pho- 
ton energy will be required to reveal the full electronic 
dispersions in the HO phase. 



D. The Fermi surface of URu2Si2 

1. Nesting vectors 

An appropriate description of the Fermi surface topol- 
ogy of URu2Si2 is an indispensable ingredient for unrav- 
eling the nature of the HO phase a s well as the uncon- 
ventional superconductivity (SC).'^™^^'^^ Experimental 
information regarding the FS of URu2Si2 has been gained 
from nesting vectors, identifled through inelastic neutron 
experiments,-® "^^ ^^ ^^^ and through extremal FS orbits , 
obtained from quantum-oscillation experiments.'^SEniinsI 

To start our discussion, we show a a side-view of the FS 
of PM URu2Si2, computed with the LSDA approach, in 
Fig. [Tsj The two FS sheets reveal the existence of a nest- 
ing vector with length c* = 27r/c (i.e., half the distance 
from one F point to the next nearest F point). These 
two FS sheets have a similar round curvature, favorable 
for nesting, with the exception that close to Z/2 the F- 




FIG. 15: (Color online) Side-view of the Fermi surface of 
PM URu2Si2 along the k^ or ky axis of the BCT BZ (see Fig. 
[I]). The arrows indicate the AF nesting vector (with length 
2tt/c) connecting the two FS sheets. Note that two smaller 
FS sheets (one Z centered ellipsoid and one small F centered 
ellipsoid) are not seen here. 



centered sheet has a more pointed part, with only a small 
area that would not be favorable for nesting. This FS 
part corresponds to a small part at the Z point in the sim- 
ple tetragonal cell, which we believe to be insignificant .I^Sl 
The identified nesting vector fits accurately to the antifer- 
romagnetic wavevector Qaf = (0; 0, 1) of longitudinal 
spin-fluctuations observed in th e HO ph ase with inelastic 
neutron scattering experiments,'2SI^^2 and it is the AF 
ordering vector of the LMAF phasc'^SEH^pj^jg nesting vec- 
tor is important for understanding the low temperature 
behavior of URu2Si2. When a coherent state emerges at 
temperatures sufficiently below the coherence tempera- 
ture T*, the system develops a FS sustaining this nesting 
vector, which is favorable for AF spin-fluctuations in the 
HO phase and formation of long-range AF order in the 
LMAF phase. Inelastic neutron experiments^ showed 
that the inelastic response at Qaf in the HO phase be- 
comes the static AF Bragg peak of the LMAF phase.1351 
The thereby induced symmetry-breaking implies a fold- 
ing of the BCT BZ at Z/2, i.e., folding Z to F. 

A second, incommensurate nesting vecto r of URu 2Si2 
has been detected at Qi = (1 ± 0.4, 0, 0).l2S132Eni This 
nesting vector has bee n o bserved in both the HO and 
LMAF phase.l^In Fig. 16 we show a cross-section of the 
LMAF and PM Fermi surfaces in the z — plane. To 
draw comparison, both FS cross-sections are plotted in 
the simple tetragonal unit cell of the LMAF phase. As 
was reported recently,'25i a clear nesting occurs (depicted 
by the dashed arrow) at 0.4 a* {a* = 27r/a), implying a 
nesting vector that matches precisely the measured in- 



commensurate vector Qi = (1 ± 0.4, 0, 0). Fig. 16 
illustrates that Qi is a suitable nesting vector for the 
LMAF phase, but less so for the PM phase (shown by 
blue circles), because in the latter phase the FS cur- 
vature does not support nesting as much. Hence, the 
incommensurate nesting vector is characteristic for the 
LMAF phase, but not for the PM phase. A priori we 
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FIG. 16: (Color online) Cross-sections of the PM (circles) 
and LMAF (squares) Fermi surface in the z = plane of the 
simple tetragonal BZ. The Fermi surface portions along the 
F— M directions are gapped in the PM to LMAF transition. 
Note that other Fermi surface parts are completely unaffected; 
the LMAF Fermi surface cross-sections are on top of the PM 
ones. The thick dashed arrow indicates the nesting vector 
0.4 a*, the thin dashed arrow indicates the position of one 
of the eight degenerate band crossings precisely at the Fermi 
energy Ef- 



don't know how the FS in the HO phase looks like, be- 
cause this would require knowing the order parameter of 
the HO, but it is known that de Haas- van Alphen exper- 
iments could not detect any notable difference between 
the FSs of the HO and LMAF phases.lSZl Likewise, inelas- 
tic neutron scattering experiments^SI give the same in- 
commensurate wavevector in the HO and LMAF phases. 
From these facts we infer that the FS of URu2Si2 in the 
HO phase should be quite close to the one we have com- 
puted for the LMAF phase. 

The gap symmetry in the superconducting phase has 
been discussed recently,HSlHIl b^t the symmetry of the 
partial gap in the HO phase has not yet been studied. 
The calculations in Fig. 



16 illustrate that the HO gap 



has in the z = plane a fourfold symmetry, in which 
there exist nodal lines with {1^2 _y2 symmetry in the ST 
structure, which is equivalent to nodal dxy symmetry in 
the BCT structure. A possible additional phase factor 
in the gap structure might exist, but cannot be deduced 
from the current calculations. 



TABLE I: Comparison of experimental and ab initio calcu- 
lated extremal FS areas in URu2Si2, for -ff || c. The FS cross- 
sectional areas are labeled according to the convention of 
Refs.l90land ll03] and have been determined with de Haas- van 
Alphen (dHv A) inea surementa^ or Shubnikov-de Haas (SdH) 
measurements.'^^'^^ The theoretical FS orbits are shown in 
Fig. |17[ All FS cross-sectional areas are given in kT. 
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dHvA 
Ref.ra 



SdH 
Ref.Umi 



SdH 
Ref.H 
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1.05 
0.42 
0.19 



1.35 
0.98 
0.50 
0.25 



0.85 
0.48 
0.22 



1.37 
0.76 
0.48 
0.26 




FIG. 17: (Color online) Computed Fermi surface sheets of 
URu2Si2 in the LMAF phase. The extremal Fermi surface 
orbits for field along the z axis are indicated. Not shown is 
the small F-centered ellipsoid (7') that is inside the large F- 
centered surface in the top panel. High-symmetry points are 
indicated in the bottom panel. 



2. Quantum oscillations 

Details of the FS of URu2Si2 has been inves- 
tigated, too, through studies of quantum oscilla- 
tions in dc Haas-van Alphen ( dHvA) and Sh ubnikov- 
de Haas (SdH) measurements. I^SlMlMliIlHill] i^ ^^^_ 



lier experiments 



[T22IH241 



only a few extremal orbits 



could be det ected. T he more recent dHvA and SdH 
measurementpSSnUnS! on purer crystals have provided a 
consistent set of data for the quantum oscillations in 



URu2Si2. These experimental data are listed together 
with the calculated extremal FS cross-section areas in 
Table |lj The dHvA experiments of Ohkuni et aL^" as 
well as the SdH experiments of Jo et a/.^^ revealed three 
FS orbits, which these authors labeled a, f3, and 7. How- 
ever, it was already notecP^lthat one or more orbits were 
missing, because the cyclotron mass of the small FS or- 
bits would not be able to account for the full, enhanced 
mass that was deduced from the specific heat. Quite re- 
cently, it was shown with SdH experiments on ultraclean 
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crystal^i23| ^j^g^^ there is indeed a further branch, which 
was named the e branch. The calculated FS of URu2Si2 
is shown in Fig. [T7| where we have also indicated the the- 
oretical extremal orbits e, a, /3, and 7 for magnetic field 
along the z axis. Note that there is a small F-centered 
ellipsoid (7') which is not visible in Fig. 17 because it 
is masked by the large F-centered ellipsoid. Table |T] il- 
lustrates that the calculated extremal orbits of the four 
branches are in quite good agreement with the available 
experiments. The computed orbits pertain to the LMAF 
phase, however, consistent with the observations made 
above, we assume the same F S fo r the HO and LMAF 
phases. The newly discovered!^ e branch corresponds 
to the large F centered rugby ball, stemming from band 
107 in our calculations. The earlier reported extremal 
orbits correspond to the M-centered ellipsoids (branches 
a and 7), the rounded half-sphere (/3), and the small, F- 
centered sphere (7'). The computed frequencies of these 
extremal orbits agree well with dHvA and SdH measure- 
ments: the extremal areas of the e and (3 orbits are well 
reproduced. The a orbit is somewhat smaller than the 
experimental value and the 7 orbit is somewhat larger. 
The calculations predict, in fact, a 7 and a 7' orbit, which 
have nearly the same frequency. Also, as discussed be- 
low, the angular dependence of their cross-sectional area 
is very similar. This similarity appears to be coincidental. 
Experimentally, it would thus be difficult to distinguish 
between the 7 and 7' branches. One other difference be- 
tween the experiments and the calculations is that the 
calculations predict an additional, small FS part around 
the Z-point (see Fig. 17 Its extremal area is very small 
(about 0.07 kT), and, as mentioned previously,'23l its oc- 
currence is related to an insignificant, small FS area at 
the Z-point which occurs through band folding. The ori- 
gin of this small FS sheet is the pointed part of one FS 
sheet near ±Z/2 in the PM phase (Fig. 15 1, which causes 
an imperfect nesting of the two FS sheets and hence an 
incomplete FS gapping near the ST Z-point when being 
folded to the LMAF phase. 

We briefly mention that earlier studies of quantum 
os cilla tions in URu2Si2 were reported by B ergemann et 
aL,li2l!ohkuni et al^a-nd Keller et aL^^Ahhough the 
measurements were performed on less well-defined sam- 
ples, there are some consistencies with the newer data. 
Bergemann et al. report the observation of two orbits, 
which cross-sectional areas of 1.09 and 0.41 kT for mag- 
netic field along the c-axis. Keller et al. reported magne- 
toresistive measurements from which they obtained sev- 
eral SdH frequencies of about 1.0 and 0.2 kT. 

The crystallographic anisotropy of the extremal FS ar- 
eas has been investigated by angular dependent dHvA 
and SdH experiments. The computed angular depen- 
dence of the FS areas is shown in Fig. [18] for field direc- 
tions in the BCT cell. A comparison of the computed an- 
gular dependence with experimental data reveals a good 
overall correspondence. Ohkuni et al^^ report that the 
a, /?, and 7 branches are quite flat, i.e., corresponding to 
nearly spherical FS sheets. Shishido et al}^^ also report 
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FIG. 18: (Color online) Calculated dependence of the ex- 
tremal FS cross-sections of URu2Si2 on the magnetic field 
direction. The extremal orbits are visualized in Fig. |17| The 
band number from which the FS sheet originates is given after 
the greek label of the respective orbits. 



that the a and j3 branches are flat and due to spherical 



FS's. As Fig. 18 shows, the theoretical FS sheets are pre- 
dicted to be relatively spherical, and hence, a very simi- 
lar angular behavior is given by the calculations. More- 
over, a more detailed comparison reveals that the angu- 
lar anisotropy of the £, a, and /3 branches is also in good 
agreement with experiment. The extremal frequency of 
the a orbit is predicted to increase from the [001] to [100] 
direction, which was also found experimentally.^- The 
frequency of one part of the j3 orbit decreases slightly, 
in accordance with experiment. The four half-spheres 
in the z = plane are anisotropic with respect to the 
field direction, therefore the single P orbit for field -ff||c 
splits in two orbits for fields along [100], but no split- 
ting occurs along [110]. Such a splitting of the / 3 br anch 
has been detected recently in SdH experiments.li^^ The 
£ orbit increases in the [100] direction, which was also 



observed in recent SdH measurements 
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However, the 



SdH experiment finds a larger increase than the theory 
predicts. 

The crystallographic anisotropy calculated for the 7 or- 
bit is not in accordance with experiment; the calculation 
predicts a slight increase of the frequency for magnetic 
fields towards [100], but the dHvA experimentPS! observed 
a slight decrease. Ohkuni et a/.l^ mention, however, that 
the signal of this orbit is very weak and therefore it was 
difficult to follow the orbit under angular field variation. 
The recent SdH experimentii^s] could not determine the 
angular dependence of the 7 orbit for the same reason. 
We also note, lastly, that the difference between the 7 
and 7' extremal orbits is predicted to be larger for fields 
along the [110] direction. 

As a nearly localized uranium 5/^ configur ation has 
been suggested to be applicable to URu2Si2,'^^'^ we 
briefiy consider the FS that would correspond to such 
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FIG. 19: (Color online) The Fermi surface of URu2Si2 in the 
LMAF phase computed with a 5f^ local uranium configura- 
tion. 
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FIG. 20; (Color online) Size of the FS gap of URu2Si2 in the 
z = plane, computed as a function of the amplitude of the 
static AF uranium moment along the c axis. 



configuration. In Fig. [19] we show the FS computed for 
URu2Si2 in the LMAF phase obtained for a localized 
5/^ configuration. Obviously, this FS is much larger than 
that computed assuming delocalized 5/'s.E^In addition, 
there are quite a number of extremal orbits (more than 
ten for i^Hz). Both these aspects of the 5p localized 
FS are not in correspondence with experiments. Con- 
versely, the overall agreement between the experimental 
and theoretical FS computed with delocalized 5/'s can 
be regarded as quite good, a result which definitely sup- 
ports that the theoretically predicted LSDA/GGA FS is 
in close agreement with the experimental Fermi surface. 



E. Longitudinal spin-fluctuations and dynamical 
symmetry breaking 

1. Potential relevance of AF spin-fluctuations 

A significant difference between the HO and LMAF 
phase is the existence of an intense mode of long-lived 
longitudinal AF fluctuations at the wavevector Qaf in 
the HO ph ase, w hich has been identified as a fingerprint 
of the HO.I^SEll This AF mode was stu died in inelastic 
neutron experiments by Broholm et al^^^, Wiebe et 
al.^and Villaume et al. ^ a.nd i n x-ray scattering exper- 
iments by Bernhoeft eiaZ.'i^Mll ^he longitudinal fluctu- 
ation was estimatecpSHSl to amount to AAI^ oc 1.2 /i^, 
its characteristic time-scale is of the order of picoseconds 
or slower .1351 The potential importance of the longitudinal 
AF mode for the HO phase was recently emphasized by 
Elgazzar et al^ 

Although the existence of the intense AF mode in the 
HO phase is evident, its connection to the HO phase is 
not yet fully understood. A major question in the dis- 
cussion is whether a mode of long-lived AF fluctuations 
can induce any FS gapping. In ordinary materials, a 



mode of (normally transversal) spin excitations would be 
a very weak perturbation of the electronic structure, and 
hence it cannot produce any significant transition in the 
bulk thermodynamic or transport properties. However, 
URu2Si2 is an exceptional material, as we find that a 
substantial FS gapping of 750 K occurs through sym- 
metry breaking on an energy scale of less than 7 K. To 
investigate the influence of the inelastic mode we consider 
the fluctuating longitudinal magnetic moment on one of 
the two uranium sublattices, i.e., Mz{t) « Mq ■ cos{ujt). 
Mq is the maximal amplitude of the longitudinal mode 
which we assume to be the same as the static LMAF mo- 
ment of 0.39 fj,B- We furthermore assume for simplicity 
a dispersionless frequency oj^q) w w. This is a good ap- 
proximation near Qaf- Obviously, the resulting average 
sublattice moment, i.e., averaged over a sufficiently long 
time-scale r, AI^ ~ t^^ J AIz{t)dt, will vanish. 

Next we consider what happens to the FS hot spots 
due to the presence of the longitudinal mode. At discrete 
time steps we can treat the AF mode as a frozen static 
antiferromagnetic arrangement with a certain amplitude 
of the moment M^- At each time snap-shot, the static 
AF arrangement will affect the FS hot spots, such that 
the degeneracy of bands at the Dirac point will be lifted. 
As the degeneracy at the Dirac point is quite sensitive, al- 
ready for small static AF moments a FS gapping appears. 



In Fig. 20 we plot the size of the FS gap as a function 
of the static AF moment. The gapping appears for small 
AF moments only at the FS hot spots, but for increasing 
moments the gap widens until the hole FS part centered 
along the F— M direction has been gapped (cf. Fig. 16 1. 



The magnitude of the computed gap develops, to a rea- 
sonable approximation, as a linear function of M^ ^ Fig. 
[20] additionally shows that the maximal FS gapping shifts 
from being at the hot spots for moderate Mz amplitudes 
to being along the F— M axis for the largest amplitude. 
For the time-dependence of the gap caused by moderately 
slow spin-fluctuations it is essential to note that the gap 
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A is an even function of Af^. Using that the magnitude 
of the calculated gap is a linear function of the moment 
Mz, we can write A(i) = a|M2(i)|, with a a positive 
constant. The time-dependence of the moment gives the 
time-dependence of the gap, A(i) — Ao|cos(a;t)|, with 
Aq = A^MAF' the maximal gap occurring for the LMAF 
phase. The time-averaged, macroscopic gap Aj^^ follows 
from 



(longitudinal) dynamical spin-spin correlation function, 



- / A{t)dt 



Ao 



\cos{ujt)\dt ^ -A^ 



(1) 
Hence, with these approximations wc estimate that the 

time-averaged HO gap is 64% of the gap of the LMAF 
phase. This value is in reasonable agreement with ex- 
perimental data obtained from resistivity measurements 
that find that the HO gap is 70 to 80% of that of the 
LMAF phase.afl-^iaoiom, 



Consequently, in the special case of URu2Si2 a long- 
lived AF mode can indeed induce a substantial FS gap- 
ping. Moreover, the similarity detected experimentally 
in the FS gapping of the HO and LMAF phases can be 
explained by the presence of the AF mode in the HO 
phase and AF order in the LMAF phase. Nonetheless, in 
spite of the similar gapping, thermodynamic and trans- 
port properties, we emphasize that the HO and LMAF 
are distinct phases. In the LMAF phase the sublattice 
U moment Mq is non-zero and thus is the standard or- 
der parameter (OP) for an ordered AF material. In the 
HO phase, conversely, M^ = and can obviously not be 
a suitable order parameter for the HO. We note that in 
the here-developed theory it is the amplitude Mq{T) of 
the longitudinal AF mode that determines the magnitude 
of the time- varying gap (cf. Fig. 20). An appropriate, 
observable secondary OP for the HO that follows from 
our model could be derived from the macroscopic aver- 
age of an even function of Mz{t). As mentioned above, 
the FS gap is such an even function, A(t) c>c |M^(<)|. 
Assuming for the temperature-dependence of the mode 
Mz{t) = Mq{T) ■ cos(a;i), this gives for the temperature- 
dependence of the averaged, macroscopic gap Aho(^) c« 
Mq{T). Consequently, if the amplitude of the mode 
Mo(T) behaves as an OP, the macroscopic, time-averaged 
gap Aho should be expected to behave as a secondary 
OP. Two experimental techniques, point contact and far- 
infrared optical spectroscopy revealed that the HO gap 
does indeed approximately behave as a typica l BCS-type 
order parameter with temperature.^^SHMlsni Vgj-y recent 
STS measurements have unambiguously demonstrated 
this feature.'''® ^"^ We remark that in contrast to standard 
BCS theory the here-obtained FS gap is not symmetric 
around Ep , see Fig. [5] 
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where Szj is the z-component of the spin opera- 
tor at the U position Ri. For relativistic materials 
such as actinides, it is more appropriate to use the 
total angular momentum, i.e., {Jz^i(t)Jzj{0)). Ap- 
proximating {Jz,i{t)Jz,j{0)) ~ (^z,j(t))(Jzj(0)), i.e., as 
Mz,i{t)Mzj{0), and using for q the AF wavevector Qj^p, 
the dynamical spin-spin correlation function S{q, lj') 
would show a resonance peak at w' = w that would 
be detectable in inelastic neutron experiments. Com- 
puting the frequency-integrated area of the peak, A = 
J S{Qjip,u}')doj', it follows that A{T) oc M§{T) oc 

2 

A (T). Hence, the area of the inelastic neutron peak 
at the AF wavevector was predicted to behave as a sec- 
ondary OP. A very recent, precise study of the inelastic 
AF resonance confirms that the area of the inelastic peak 
indeed behaves as a BCS-type OP.I^In the LMAF phase 
the AF mode seizes to exist and hence the inelastic peak 
area vanishes. Adopting the spin-spin correlation func- 
tion expressed through the inelastic peak area A as a 
derived OP for the HO phase, one obtains the situation 
where the staggered AF moment M^ ^ and A = m 
the LMAF phase, while conversely Mz = and A ^ in 
the HO phase, i.e., the phases have distinct OP's (cf. the 
discussion regarding this in Refs. [T| l33l [84l I13H I132| and 



A second possibility for an even function of Mz{t) 
would be the dynamical susceptibility, expressed by the 



Interestingly, so far only the gap Ajjq and the inelastic 
neutron pea k inten sity A were proven to display typical 
OP behavior .E5M1 Both quantities are evidently related 
to the intense AF mode in the HO phase, emphasizing 
that the AF mod e is essential for the HO. Earlier ^^Si 
NMR experiment J^^ indicated that the NMR line width 
displayed OP behavior in the HO phase. Another NMR 
experiment, however, observe d an inhomogeneous devel- 
opment of AF order below Tq ^^ and a further NMR ex- 
periment did indeed detect an increase of the line width 
below Tq, but this effect disappeared at lower tempera- 
tures and more strain-free crystals.!^ 

Second order phase transitions are characterized by a 
divergence of a (generalized) susceptibility at the crit- 
ical temperature. We note that the critical behavior 
of the dynamical susceptibility x{Qt^) ^t the HO tran- 
sition is consistent with the proposed importance of 
spin-fiuctuations. Its imaginary part x" is proportional 
to the spin-spin correlation function, S{q,uj). Its real 
part x'(q,w), however, was found-from inelastic neutron 
experiments-to exhibit a s harp cusp at q = Qaf and 
a; w when entering the HO.^i^This signals a divergence 
of x'(Qaf, w ~ 0) at To, which has been broadened only 
by experimental resolution. 

The magnetic fluctuations could also bear relevance^Sl 
to the occurrence of unconventionaP^"^ SC out of the 
HO at temperatures below 1.2 K.EHSl The persistence 
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of spin-fluctuations at the incommensurate wavevector 
d own i nto the SC phase was observed by Br oholm et 
q;j 39 | 42 | ^ very recent inelastic neutron studjff^ reveals a 
small energy shift in the AF mode upon entering the SC 
phase, suggesting an involvement of the AF spin excita- 
tions in mediating the SC. 



2. Symmetry breaking 

Static AF ordering obviously corresponds to break- 
ing of time-reversal symmetry and of the body-centered 
translation vector in the LMAF phase; the latter causes 
a doubling of the unit cell along the c axis. What would 
be the symmetry breaking induced by the longitudinal 
AF mode? The unconventional answer is, that it in fact 
depends on the time-scale being considered. On a very 
short time-scale, typically much less than that of one os- 
cillation, the body-centering as well as time-invariance 
will "instantaneously" be broken. On a sufhciently long 
time-scale however the moments on both uranium atoms 
have reverted their orientation many times, and thus the 
two atoms are no longer distinguishable. As the two 
U atoms are undistinguishable on long time-scales there 
will, consequently, be no breaking of time-reversal sym- 
metry nor of the c-axis translation vector. Alternatively, 
this can be assessed, for example, from the moment- 
moment correlation function (Eq. (|2|), which is time- 
even. Most experimental techniques probe quantities on 
relatively long time-scales, much longer than picoseconds 
(e.g., ^SR or elastic neutrons). These will not be able to 
detect a symmetry breaking in the HO phase. 

The question whether time-reversal or possibly trans- 
lational symmetry is broken in the HO is an issue of 
on-going debate. Earlier neutron experiments provided 
evidence for the breaking of time-reversal symmetry,'^ 
however these measurements might have been influenced 
by the presence of parasitic small moments. Also later, 
more precise neutron experiments deduced evidence for 
time-invariance breaking,ESiIi32l yet the issue is not un- 
ambiguously settled. Newer ^^Si NMR on single crystals 
detected a small line shift as well as a linewidth enhance- 
ment below Tq which corres pond s to the presence of a 
small internal magnetic fleld.'^^SI This flnding might in- 
dicate the presence of time-invariance breaking in the 
HO phase. Also, the temperature-dependent static mag- 
ne tic suscep tibility of URu2Si2 displays a clear kink at 
j|^j ii | i39 [i40] gmjij feature is not uncommon at an AF or- 
dering transition in actinide compounds, and it indicates 
the influence of a magnetic field on the HO. This in turn 
could suggest that a purely electric OP such as an electric 
quadrupole or hexadecapole is unlikely. Electric hexade- 
capolar order has been proposed for PrRu4Pi2 precisely 
because of the observed i nsen sitivity of an ordering tran- 
sition to a magnetic fleld.QSll 

Altogether, clear experimental evidence for cither 
breaking of the time-reversal or the translational sym- 
metry in the HO phase is so far lacking. A possible 



verification of the symmetry-breaking due to the AF 
mode would require a sufficiently fast experimental tech- 
nique, such as x-ray scattering or photoemission. The 
here-proposed breaking of the body-centering could be 
probed in ARPES experiments, where the correspond- 
ing doubling of the unit cell would appear as a sudden 
folding of bands at Tg, together with a change in the c* 
axis periodicity. A further argument in favor of time- 
reversal and translational symmetry breaking in the HO 
might come from dHvA and SdH experiments. The the- 
oretical FS (Fig. 17) computed for the symmetry-broken 
phase agrees well with the inferred experimental FS. This 
correspondence suggests that the body-centered transla- 
tion vector is broken in the HO phase and most likely 
the time-reversal symmetry as well. Conversely, using 
the FS computed for the bct phase, it is not possi- 
ble to exp lain the observed dHvA and SdH quantum 
oscillations I 36 | 90 | i 03| Although this is consistent with the 
dynamical symmetry breaking model, it cannot be ex- 
cluded that other mechanisms for the HO lead to the 
same symmetry reduction. 



V. DISCUSSION 

It is elucidative to consider if hidden order phases 
have been discovered in other 5/ materials. Apart from 
URu2Si2, the actinide oxide Np02 has drawn consider- 
able attention,EI22Hilll because a sharp phase transition 
to an unusual phase occurs below 30 K as is witnessed 
by specific heat measurements.^"*^ Intensive research dur- 
ing the last two decades revealed that, in the absence of 
any magnetic dipole moment, long-range multipolar or- 
dering of a higher-order magnetic multipole (octupole or 
higher) on the N p ion s is likely to occur in the unusual 
phase below 30 K.EISSlxhese higher multipoles cannot be 
observed experimentally, but electric quadrupolar order 
appearing as a secondary OP has been observed through 
resonant x-ray scattering (RXS)"^ and ^^O NMR.^^ 
Theories based on 5/'^ localized states on the Np**"* " ion 
have further provided insight in the multipolar order.'^'^^ 

Also for URu2Si2 a number of theories for the low- 
temperature order have been developed on the basis of 
localized or nearly localized 5/^ configurations. Specif- 
ically, for URu2Si2 the following mul tipolar OPs have 
been pr oposed: electric quadrupole,l2^ElIlo| magnetic 
octupole,'i2l21I23lIi] electric hexadecapole,'2£l and magnetic 
triakontadipole.'23 Other theories on the basis of crystal 
electrical field (CEF) consider ations for a U'*"'" 5/^ ionic 
state have also been developed.C^Ell Electric quadrupolar 



ordering should be detectable withRXS,'^^^ yet experi- 
mental studies gave a null result .'^^^'i^ Quadrupolar or- 
der requires an Ei — Ei scattering channel to be detected, 
whereas higher order multipoles require at least an E2 
optical transition. Since the scattering cross-section for 
E2 transitions is extremely small, a detection of higher 
order multipoles is currently unlikely. Furthermore, most 
of these theories are based on a localized U 5/^ configu- 
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ration, but in spite of a similarity in the thermodynamic 
properties of URu2Si2 and Np02, treating the 5/'s as 
locaHzed in Np02 is a good approximationfi^ but it is 
questionable if this approximation is valid for URu2Si2, 
too. 

Triakontadipolar order in the HO phase of URu2Si2 
was recently proposecP^ on the basis of LSDA+[/ cal- 
culations with a Coulomb U of about 1 eV. In con- 
trast to localized 5/^ theories, the / electrons are rel- 
atively delocalized in this description. These LSDA+C/ 
calculations predict an unusual long-range ordered AF 
state for URu2Si2 with parallel s]iin and orbital magnetic 
moments. ^'^^ URu2Si2 is an unconventional material and 
it might indeed be that such an unusual magnetic state 
is realized. It should be noted however that thus far not 
a single magnetically ordered U-compound has been dis- 
covered having this property. This can be understood 
from the strong spin-orbit interaction in actinides, which 
enforces always antiparallel spin and orbital moments 
in the early actinides. Also, conventional LSDA calcu- 
lations predict antiparallel spin and orbital moments,^* 
and, as mentioned above, antiparallel spin and orbital 
moments are in fact consistent with neutron form factors 
measurements. 'SS Moreover, the values deduced for the 
separate spin and orbital moments are in close agreement 
with values predicted by LSDA calculations. The Fermi 
surface predicted by LSDA-I-C/ calculations deviates al- 
ready substantial from the LSDA FS. In particular, new 
band crossings appear and the FS gapping feature is lost 
(cf. Fig. 10); thus, the LSDA-fC/ approach with a f7 of 



more than 1 eV cannot explain the FS gapping appeari ng 
in the LMAF phase. The recent LSDA+t/ calculationdSl 
predict that, with expansion of the lattice parameter a, a 
phase transition would occur in URu2Si2 from the AF or- 
dered phase to the same AF ordered phase, but with both 
spin and orbital moments being opposite to the original 
ones. The triakontadipolar moment is predictecP^l to be 
present in both AF phases, but it could become the main 
OP precisely at the transition point between these two 
AF phases, where the dipolar moment would vanish. It 
should be noted, however, that such a phase transition 
from one AF to another AF phase, being identical by 
symmetry has not been observed in URu2Si2. Also the 
predicted^^ ^^^ parallel spin and orbital moments are not 
supported by experimental observations. 

Electric hexadecapolar order in URu2Si2 has recently 
been proposed on the basis of DMFT calculations.^^ 
These calculations suggest a nearly localized 5/^ con- 
figuration on the uranium ion, consisting of two singlet 
CEF levels, separated by 35 K. The HO phase emerges 
from an excitonic mixing of the CEF ground and the 
first excited state singlet. The same CEF level scheme 
(but notably with gro und and excited state reversed) 
had been proposecP^ earlier to explain the Ising-like 
anisotropy and maximum of the magnetic susceptibility. 
Although such CEF model can explain certain properties 
of URu2Si2, it is one of the observations of the present 
study that the CEF model does not corroborate with 



many other experimental data. 

The here-reported computed results suggest alto- 
gether that an applicable explanation of the prop- 
erties of URu2Si2 ought to arise from an itinerant 
5/ ele ctronic structu re. A number of earlier HO 
theorie^ ^ ' ^''' ' ^^ ' ^" ' ^^ -^^ have been based on the itinerant 
5/ picture, and more recently, two theories have empha- 
size d the importance of dyna mic spin-fluct uations for the 
jjO.mmi The earher theoriePE^HlSoIlll ^^^^A offer an 

explanation for some aspects of the HO, but none of 
the models could be unambiguously confirmed experi- 
mentally (see, e.g., Refs. |3H1 IZHl and 11541 for a discus- 
sion). The two HO theories that highlight the role of 
dynamical spin excitations have either focussed on spin- 
fluctuations at the AF wavevector^S or at the incom- 
mensurate wavevector.l^ Obviously, spin-fluctuations at 
the incommensurate wavevector as proposed in the re- 
cent theory by Balatsky et alW^ are expected to have 
an influence, therefore further investigations of the in- 
commensurate mode are needed to establish the relative 
importance of the two modes. 

In the present work we have studied particularly the ef- 
fect of the AF mode and find that its presence in the HO 
phase can explain the FS gapping, the broken-symmetry 
FS, the relation A^^ sa (0.7 — 0.8)Aj^„^p, as well as 
the mean-field OP behavior of A^^^ and of the inelastic 
neutron peak. The dynamical symmetry breaking model 
nonetheless builds on the surprising and exceptional ex- 
planation of the HO phase being driven by the dynamic 
AF mode. Spin excitations are usually only weak pertur- 
bations of the ground state, that in conventional mate- 
rials, cannot modify the ground state nor its properties. 
In URu2Si2, however, the exceptional situation appears 
to be realized that the low-lying spin excitations actually 
dictate the thermodynamic properties of the HO phase. 
Thus far a similar situation was app arently realized only 
in one other material, PrAu2Si2,'^^ which has the same 
crystal structure as URu2Si2. Why the low- lying AF spin 
excitations are so effective in URu2Si2 is related to mate- 
rials' specific aspects. The calculated energy scale of AF 
excitations is only of the order of 7 K, but they couple 
to an unexpectedly large FS reconstruction, with gaps of 
about 700 K. Thereby these low-lying modes can essen- 
tially modify the thermodynamic and transport proper- 
ties and are indeed inherent to the HO phase. 

Our calculations lastly emphasize that in URu2Si2 
several remarkable, materials' specific features are com- 
bined. First, the LMAF phase with its relatively mod- 
est total moment has a total energy that is very close 
to that of the PM phase. Second, the low-temperature 
FS of URu2Si2 supports a nesting vector Qaf that pro- 
motes long-range AF ordering or longitudinal AF spin- 
fluctuations. Third, the energy bands in the ST BZ ex- 
hibit accidental degenerate crossing points at Ep as well 
as close to Ep. Fourth, a breaking of the c axis trans- 
lational symmetry together with time-reversal symme- 
try breaking in the LMAF phase or through longitudinal 
spin-fluctuations causes a lifting of the degeneracy at the 
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Dirac points and thus to the opening of a partial gap. 
These observations underhne that to a large extent the 
HO is not a generic property, but rather, it is borne out 
of the materials' specific electronic structure. 



VI. CONCLUSIONS 

For a quarter century the mysterious hidden order 
phase in URu2Si2 has been in the focus of many in- 
vestigations. A detailed understanding of the underly- 
ing electronic structure of URu2Si2 in both the PM and 
LMAF phases is required before a full explanation of the 
HO can be formulated. Neither the PM nor the LMAF 
phase appear to be exceptional, therefore first-principles 
electronic structure calculations should be capable of ex- 
plaining the known solid-state properties of these two 
phases. Using ab initio electronic structure calculations 
on the basis of the DFT-LSDA/GGA, the GGA+U, and 
LDA+DMFT methodologies, we have performed an ex- 
tensive study of the electronic structure of URu2Si2. A 
major question has for a long-time been what the appli- 
cable electronic structure of the PM and LMAF phases is, 
whether it has itinerant or localized uranium 5/ electrons 
or perhaps 5/ electrons with dual, i.e., both itinerant and 
localized characteristics (cf. Ref. I156p . Our conclusion 
regarding this issue is that we obtain an electronic struc- 
ture picture consistent with available experimental data 
only when we adopt delocalized 5/ states. 

Specifically, we find that materials' specific DFT- 
LSDA/GGA electronic structure calculations on the 
basis of delocalized 5/'s explains the following low- 
temperature properties of URu2Si2: 1) the equilibrium 
volume, 2) the internal zsi coordinate and c/a ratio, 3) 
the bulk modulus and equation of state, 4) the spin and 
orbital magnetic moment of the LMAF phase, 5) the 
closeness in total energy of the PM and LMAF phases, 
with the LMAF phase becoming more stable under pres- 
sure, 6) the Fermi surface gapping and instability, 7) the 
compensated metal character, 8) the number of holes, 9) 
the 5/ occupation number, 10) the resistivity jump at the 
PM to LMAF transition, 11) crystallographic anisotropy 
of the resistivity, 12) the gapping in the infrared optical 
spectra, 13) the antiferromagnetic and incommensurate 
FS nesting vectors, and 14) the dHvA and SdH frequen- 
cies, as well as their angular dependence. In contrast, 
assuming a localized uranium 5/^ configuration we can- 
not obtain a satisfactory explanation of the experimental 
low-temperature data. The recent ARPES datgt^' form 
a notable exception to the set of low-temperature prop- 
erties that are explained by 5/ itinerant electronic struc- 
ture calculations. A (nearly) localized 5/^ uranium con- 
figuration appears to tally better with the observed spec- 



tral features.!^ As the Fermi surface detected through 
ARPES evidently has to match that obtained with quan- 
tum oscillation techniques, more experimental investiga- 
tions appear to be needed to resolve this issue. 

We have additionally studied URu2Si2 in the high- 
temperature PM phase, using LDA-(-DMFT calculations. 
Our DMFT calculations predict the progressive opening 
of a quasi-particle coherence gap at the chemical poten- 
tial when temperature is reduced toward the coherence 
temperature. 

An explanation of the hidden order requires still taking 
a step beyond conventional electronic structure calcula- 
tions, however, the electronic structure picture underly- 
ing any explanation of the the HO phase must be consis- 
tent with those underlying the PM and LMAF phases. 
The current investigations strongly emphasize that a 
well-grounded explanation of the HO has to emerge from 
an itinerant /-electron picture. From the presented cal- 
culations we conclude that the presence of the intense 
AF mode in the HO phase could explain (i) the FS gap- 
ping occurring in the HO phase, (m) its broken-symmetry 
FS extremal orbits and incommensurate nesting vector, 
(Hi) the experimentally found relation between the gaps 
in the HO and LMAF phases, A^^ ss (0.7 - 0.8) A^^^^^, 
and (iv) the OP behavior of the gap Ajj^ and the in- 
elastic dynamical susceptibility. The dynamical symme- 
try breaking model appears so far to be the only theory 
that explains and even predictecP^ these four properties. 
The importance of this mode for bringing about the HO 
transition is further exemplified by its occurrence and 
gapping in its magnon dispersion, its assistance in the 
entropy removal, as well as the divergent behavior of the 
staggered susceptibility x' at Tq . The presence of the in- 
elastic AF mode at low-temperatures, down into the SC 
phase,'i2II could provide a clue as to why an unconven- 
tional form of superconductivity develops. 
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